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ABSTRACT 

The mitochondrial genome exists in numerous 
structural conformations, complicating the study 
of mitochondrial DNA (mtDNA) metabolism. Here, 
we describe the development of 2D intact mtDNA 
agarose gel electrophoresis (2D-IMAGE) for the sep- 
aration and detection of approximately two-dozen 
distinct topoisomers. Although the major topoiso- 
mers were well conserved across many cell and 
tissue types, unique differences in certain cells 
and tissues were also observed. RNase treatment 
revealed that partially hybridized RNAs associated 
primarily with covalently closed circular DNA, con- 
sistent with this structure being the template for 
transcription. Circular structures composed of 
RNA:DNA hybrids contained only heavy-strand 
DNA sequences, implicating them as lagging-strand 
replication intermediates. During recovery from 
replicative arrest, 2D-IMAGE showed changes in 
both template selection and replication products. 
These studies suggest that discrete topoisomers 
are associated with specific mtDNA-directed 
processes. Because of the increased resolution, 
2D-IMAGE has the potential to identify novel 
mtDNA intermediates involved in replication or tran- 
scription, or pathology including oxidative lineariza- 
tion, deletions or depletion. 

INTRODUCTION 

The mammalian mitochondrial genome is present in 
hundreds to thousands of copies per cell and encodes es- 
sential components of oxidative phosphorylation 
(OXPHOS) complexes I, III, IV and V. Mitochondrial 



DNA (mtDNA) is located within the mitochondrial 
matrix in single or multi-genome protein-DNA structures 
termed nucleoids (1-6), often associated with the inner 
membrane (7,8). Total mtDNA copy number is under 
tight cellular control and responds to changes in demand 
for mitochondrial respiration [reviewed in (9)]. The main- 
tenance and integrity of this genome is vital, as alterations 
such as point mutations, insertions and deletions, or low 
copy number (depletion) impair mitochondrial respiration 
and result in adenosine triphosphate insufficiency that 
can cause a wide spectrum of heterogeneous multi- 
systemic disorders (10-13). Consistent with their causative 
role in disease, mtDNA abnormalities are increasingly 
associated with a range of human pathologies, including 
neurodegenerative disorders (14,15), type 2 diabetes 
mellitus (16-18), cancer (19-21) and aging (22,23). 

Mitochondrial function and mtDNA integrity are main- 
tained by a dynamic balance among replication, transcrip- 
tion, repair and degradation. It is not yet known how this 
balance is regulated, or how specific mtDNA nucleoids (or 
genomes within the nucleoid) are assigned to each process. 
It is possible that genomes performing different roles 
may be structurally distinct; however, the inability to 
isolate mtDNA structures has limited the study of these 
processes. We expected that the separation and detection 
of mtDNA based on size, shape, double or single-stranded 
content and topological conformation would identify mol- 
ecules involved in specific aspects of mtDNA metabolism. 
To begin to address this hypothesis, we first explored 
established methods to separate mtDNA structures 
based on their shape. 

Agarose gel electrophoresis is a robust method for 
separating plasmid DNA topoisomers, and it was first 
used to examine mtDNA damage in isolated HeLa cell 
mitochondria (24,25). Improved ID separation methods 
can detect changes in mtDNA topology caused by transi- 
ent knockdown or induced expression of the mtDNA 
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organizing protein and mitochondrial transcription factor 
A (TFAM) in mouse cells (26), and they were later used to 
identify complex tangled structures (27), catenanes and 
double-length linear molecules (28) present in human 
heart mitochondria. The addition of a second gel separ- 
ation step, which uses ethidium bromide (EtBr) intercal- 
ation to compact and accelerate topologically confined 
DNA, can further resolve undigested genomes, but it 
has only been applied to petite yeast mtDNA (29) and 
high-molecular weight structures in heart (27). This 
approach contrasts the classic 2D agarose electrophoresis 
gels (2D-AGE) used to examine fragments of genomes 
(30). In the present study, we adapted the concept of 
resolving undigested (intact) mtDNA to optimize the sep- 
aration of all major structures, a technique that we call 2D 
intact mitochondrial DNA agarose gel electrophoresis 
(2D-IMAGE). We have used a detailed molecular 
biology approach to describe and characterize the struc- 
tures revealed by 2D-IMAGE, many of which are 
associated with replication and possibly transcription. 



MATERIALS AND METHODS 

Cell culture and treatments 

Mouse embryonic fibroblasts (MEFs), C2C12 myoblasts, 
MIN6, HeLa, HMEC 151p and MDA-MB-231 cells were 
cultured in Dulbecco's modified Eagle's medium supple- 
mented with 25 mM of glucose, 1 mM of sodium pyruvate, 
lx pen/strep and 10% of fetal bovine serum:fetal calf 
serum (1:1). MIN6 cell medium was further supplemented 
with 0.001% of 2-mercaptoethanol. SKBR3 cells were 
cultured in Roswell Park Memorial Institute medium sup- 
plemented with 10% of fetal bovine serum. Isolation of 
cardiomyocytes from neonatal myocardium was per- 
formed by the myocyte core of the Penn Cardiovascular 
Institute at the University of Pennsylvania using standard 
procedures (31). 

EtBr treatment of wild-type MEFs was conducted on 
subconfluent cells for 24 h in the presence of 1 mM of 
pyruvate, 0.2 mM of uridine and either 100, 250 or 
500ng/ml EtBr. For EtBr recovery experiments, 
subconfluent MEFs were grown in 500ng/ml EtBr plus 
pyruvate and uridine for 24 h, at which time the medium 
was aspirated, cells were washed once with phosphate 
buffered saline without magnesium or calcium (MP 
Biomedicals) and fresh medium lacking EtBr was added. 
Cells were harvested at 1, 3 and 6h after medium change. 
Cell collection was performed as follows: cells were 
trypsinized, neutralized with media plus serum, pelleted 
at 200 RCF for 5min, washed with phosphate buffered 
saline, collected at 4°C by centrifugation, flash frozen in 
liquid nitrogen, then stored at — 80°C for later DNA 
preparation. 

Relative mtDNA content by quantitative polymerase 
chain reaction (qPCR) was conducted using MT9 and 
MT11 primers (32) for mtDNA and NDUFV1 F and 
NDUFV1 R primers for nuclear DNA (33) in a 
RealPlex PCR mastercycler (Eppendorf) with SYBR 
green mastermix (Applied Biosy stems). 



DNA isolation, Southern blotting, 1-D and 2-D 
electrophoresis 

Please see Supplementary Methods for description of en- 
richment of mitochondria, DNA preparation, Southern 
blotting, ID and 2D electrophoresis and step-by-step 
protocols. 

DNA treatment 

For enzymatic treatment of DNA, samples were diluted 
into the appropriate buffer provided by each manufac- 
turer and treated with the indicated units of enzyme at 
the following temperatures and times: topoisomerase I 
(Invitrogen; 20 U at 37°C for 1 h); topoisomerase II 
(USB Affymetrix; 40 U at 37°C for 1 h); topoisomerase 
IV (Inspiralis; 20 U at 37°C for 1 h); DNA gyrase (New 
England Biolabs; 10 U at 37°C for 1 h); Bglll (Fermentas; 
4U at 37°C for 30min); SI nuclease (Promega; 0.9 U at 
37°C for 30min); Escherichia coli exonuclease I (New 
England Biolabs; 20 U at 37°C for 30min); RNase A 
(Ambion; 2ul at 24° C for 1 h, used to generate data in 
Figure 4 and Supplementary Figures SI and S2); RNase 
H (Fermentas; 15 U at 37°C for 30min). 

RESULTS 

ID separation of mtDNA reveals major topoisomers 

We first examined the qualitative distribution of total 
mtDNA topoisomers in wild-type MEFs or C2C12 
myoblasts by an established ID agarose gel-based 
method (26). RNase A-treated total wild-type MEF 
DNA (10 ug based on OD 260 ) was electrophoresed 
through a low-percentage agarose gel overnight at low 
voltage to allow higher molecular weight DNAs to 
resolve (see Supplementary Methods for details). DNAs 
were transferred to nylon membranes and Southern 
blotted for mtDNA with a probe for mouse COXI 
sequence (Supplementary Figure SI). Catenated circles, 
nicked circles, linear molecules and supercoiled circles 
were resolved in a ladder pattern (Figure 1A), typical for 
this approach (25,26,28,34-37). Preliminary characteriza- 
tions to confirm the predicted identities of the major bands 
in each lane included treatment of total MEF DNA 
with the Bglll restriction enzyme, which cleaves mouse 
mtDNA once, and topoisomerase I, which relaxes super- 
coiled DNA but does not separate catenated molecules 
(Figure 1A). DNA molecules sensitive to linearization by 
Bglll treatment are presumed to have a circular form and 
to be double-stranded (ds) at the restriction site; therefore, 
the location of the Bglll linearization product marks 
the position of genome-length linear mtDNA in the 
undigested ID gel profile. Several other low abundance 
molecules were detected after restriction digest (each 
indicated by asterisks) and were addressed in later experi- 
ments (see Figure 2). The position of maximally super- 
coiled mtDNA molecules was revealed by topoisomerase 
I treatment, which converted them to relaxed circles 
(Figure 1A and B). To confirm the location of supercoiled 
molecules, we also treated samples with DNA gyrase, 
which increases the degree of supercoiling of closed 
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Figure 1. Separation of mtDNA into major topoisomers: catenanes, relaxed circles, linear molecules and supercoiled circles. (A) Treatment of total 
cellular DNA isolated from MEFs with Bglll and topoisomerase I (Topo I). Untreated control (untr.) is shown for reference. Additional structures 
present after Bglll digest are indicated by asterisks. (B) Treatment of total cellular DNA isolated from C2C12 myoblasts with topoisomerases I, II, 
IV and gyrase. (C) Treatment of total cellular DNA isolated from C2C12 myoblasts with SI nuclease and E. coli exonuclease I. All gels were probed 
by Southern blot with COXI mtDNA sequence. 
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Figure 2. 2D separation of mtDNA resolves previously undetected structural forms. (A-F) Total MEF DNA resolved by 2D-IMAGE and detected 
by Southern blot analysis with total mtDNA probe: (A) untreated; (B) DNA from isolated mitochondria; (C) topoisomerase I treated; (D) topo- 
isomerase II treated; (E) Bglll digested; and (F) SI nuclease treated. Within panel (A) above 2D-IMAGE profile is a ID lane in the correct 
orientation for the second dimension separation for comparison, with the following annotated species: catenanes (cat.); unit length nicked circles 
(In NC); two-genome length nicked circle catenane (2n NC); CCCs; unit-length linear molecules (In L). Arrow in panel A indicates position of CCC 
bar, which collapses by topoisomerase I or II treatment (black triangle in panels C and D). In panel D, the signal reduction in catenanes is 
completely recovered in linear molecules, all signal from CCC bar collapses onto relaxed CCC and a small amount of mtDNA at the 2n NC 
region that is resistant to topoisomerase II treatment is indicated by number sign. In panel E, single asterisks indicate structures potentially also 
identified in Figure 1 after Bglll digestion, whereas open diamonds indicate additional molecules that were undetected in Figure 1. In panel (F), 
single-stranded species removed by SI nuclease treatment are indicated by black squares. 
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circular DNA (Figure IB). DNA gyrase treatment did not 
promote supercoiling of all relaxed circles, suggesting that 
the relaxed circles were not exclusively a covalently closed 
population. We next treated total genomic DNA with type 
II topoisomerases (topoisomerase II and IV), which 
decatenate or relax supercoiled molecules but leave 
concatemers (head-to-tail dimers) intact. Both type II 
topoisomerases relaxed supercoiled molecules and 
resolved catenated circles (Figure IB). 

In addition to the four major topoisomer species 
resolved by ID electrophoresis, we observed significant 
hybridization to species migrating between each of the 
major bands within each lane. Some of this signal corres- 
ponded to low-level species that were not cleaved by re- 
striction digestion and, therefore, could contain 
single-stranded (ss) DNA (Figure 1A). The less 
abundant species, including multiple bands between the 
relaxed and supercoiled circles, between relaxed circles 
and 2n catenanes, and above 2n catenanes, have been 
observed by others (26,34). To test whether any of these 
structures contained ssDNA, we treated total DNA with 
two common enzymes that remove unpaired DNA 
(Figure 1C). Treatment with low levels of SI nuclease, 
which digests ssDNA, dsDNA at nicks or gaps and 
RNA [with 5-fold lower activity than for ssDNA (38)], 
reduced the signal that was not associated with the 
major topoisomer s and produced more relaxed circles 
and various sizes of linear mtDNA. Treatment with 
E. coli exonuclease I produced no clear changes in the 
resolved molecules, suggesting that the ssDNA did not 
contain 3 / -overhangs. Although the overall pattern in ID 
gels agreed with previous findings that mtDNA is present 
in multiple species (24-26,39), we observed numerous 
other hybridizing species, including those around relaxed 
circular, linear and supercoiled mtDNA molecules. These 
data show that ID gel electrophoresis is insufficient to 
separate and identify all mtDNA topological and struc- 
tural isomers, including those that apparently contain 
ssDNA. 

2D separation of mtDNA reveals ~25 structures, 
including ssDNA molecules 

To resolve and characterize these additional mtDNA 
species, we developed 2D-IMAGE. This procedure intro- 
duces a second electrophoretic separation step after the 
ID gel that uses EtBr to induce supercoiling in covalently 
closed or topologically constrained DNAs, resulting in 
their accelerated migration. Furthermore, EtBr enhances 
the rigidity of DNA, thus exaggerating separation of mol- 
ecules with differences in extended DNA shape. After elec- 
trophoresis in the first dimension, the gel was incubated in 
EtBr, slices were excised and rotated 90°, embedded in a 
second gel of the same percentage that contained EtBr and 
samples electrophoresed under the same conditions as the 
first dimension (step-by-step method is detailed in 
Supplementary Methods). Because the agarose concentra- 
tion is constant in both dimensions, distinct populations 
of DNA appear in the shape of the well used in the first 
dimension. The four major topoisomers resolved in ID 
gels (catenanes, relaxed circles, linear molecules and 



maximally supercoiled circles) migrated near the 
diagonal of the 2D gels, indicating that their mobilities 
were largely unaffected by EtBr; however, numerous add- 
itional species were affected by EtBr and were now 
separated from these major structures. DNA from 
enriched mitochondria produced a similar pattern, 
eliminating the possibility that these structures are hybrid- 
ization artifacts of using whole-cell DNA preparations 
(Figure 2 A and B). 

Among the new DNA structures visible in the 
2D-IMAGE profile was a horizontal continuum of 
signal spanning the region between supercoiled and 
relaxed circular DNA in the first dimension (Figure 2A, 
arrow). These molecules are covalently closed circles 
(CCCs), which are separated in the first dimension by 
writhe, but have equal migration characteristics in the 
second dimension because of the presence of EtBr, 
which drives closed circles into a maximally supercoiled 
state. These molecules ranged from completely relaxed 
[left-most molecules that co-migrated with open (i.e. 
nicked) circles in the first dimension] to maximally super- 
coiled (right-most molecules corresponding to the species 
labelled 'supercoiled' in the first dimension). The CCCs 
are clearly separated from unit-length nicked circular 
DNAs (In NC), whose migration is EtBr insensitive. 
The identity of CCCs was confirmed by topoisomerase I 
treatment (Figure 2C), which converted all supercoiled 
CCCs (scCCC) to relaxed CCCs (rCCC) that still 
remained EtBr sensitive in the second dimension 
(indicated with black triangle). These data reveal a 
broad population of supercoiled molecules not detected 
by ID methods. 

To remove catenanes, we treated total genomic DNA 
with topoisomerase II, which after separation by 
2D-IMAGE revealed a small number of remaining mol- 
ecules that are likely concatemers (Figure 2D, indicated by 
number sign). The data suggest that the topoisomerase 
reactions were complete because scCCC was fully con- 
verted to rCCC. The increase in linear forms was initially 
unexpected, but has been observed previously for topo- 
isomerase II (40). 

Linear mtDNAs (In L) were identified by digestion with 
Bglll to collapse all full-length circular dsDNA 
(Figure 2E). A diagonal trail of signal originating at the 
linear molecules and proceeding towards lower molecular 
weights represents linear mtDNAs of decreasing length 
(these are included in quantitations of all linear molecules 
in subsequent experiments; see Supplementary Figure S2). 
As observed in ID analyses, restriction digestion con- 
verted the majority of mtDNA to linear molecules of 
unit length, yet several species remained. Those structures 
detected in ID (Figure 1A) are labelled with an asterisk in 
Figure 2E, whereas those only detected by 2D-IMAGE 
are indicated with an open diamond. 

To test whether any of the additional species present 
after Bglll digest contained ssDNA, we treated wild-type 
MEF DNA with SI nuclease (Figure 2F) and observed a 
complete disappearance of the two near-vertical series of 
molecules proceeding downwards from In NC and linear 
molecules and of the haze emanating diagonally from the 
linear molecules towards the origin, all of which are 
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Figure 3. Model for mtDNA structures resolved by 2D-IMAGE. (A) 2D-IMAGE profile of DNA isolated from MEF mitochondria in Figure 2B. 
Each population of mtDNA is indicated with a number. (B) Diagram of molecules resolved by 2D-IMAGE. The numbers correspond to populations 
indicated in panel (A). Those numbers with annotations are predicted molecular species based on enzyme sensitivity and electrophoretic character- 
istics. The four primary topoisomers are indicated in bold. Species 26 is found in other cell types and tissues but not MEFs. Positions of nicks in 
DNA are arbitrary. Abbreviations include: catenane (cat.); nicked circle (NC); concatemer (CNC); intermediate between NC and cat. (ln-2n); linear 
(L); covalently closed circles (CCC); relaxed CCC (rCCC); supercoiled CCC (scCCC); circle (C); single-stranded (ss). Lengths are indicated relative to 
one genome unit length (In). 



Table 1. Cell lines and tissues used in this study 



Mouse cell lines 




MEF 


Embryonic fibroblast 


MIN6 


Pancreatic insulinoma 


C2C12 


Myoblast 


MC3T3-E1 


Pre-osteoblast 


Human cell lines 




HMEC 151P 


Normal mammary epithelium (hTERT) 


MDA-MB-231 


Mammary epithelial cancer 


HeLa 


Cervical cancer 


SKBR3 


Mammary epithelial cancer 


Mouse tissues 




Liver 




Heart 




Quadriceps 




Kidney 




Human tissues 




Pancreatic islets 




Myocardium 





marked with black squares. We also detected some 
'sharpening' of the molecules between In NC and 2n 
NC, which we term ln-2n molecules. The SI -sensitive 
species were not affected by E. coli exonuclease I treat- 
ment (data not shown). Taken together, these results 
identify 25 mtDNA species (summarized in Figure 3) 
and suggest that 2D-IMAGE could be applied to the 
study of the biogenesis and metabolism of these molecules. 



Similar mtDNA structural forms exist across human and 
mouse cell and tissue types 

The diversity in mtDNA structures separated by 
2D-IMAGE analysis prompted us to determine how 
well-conserved the mtDNA structural profiles might be 
in several human and mouse cell lines. We isolated total 
DNA (treated with RNase A) from four lines each of 
mouse and human cells listed in Table 1 and subjected 
10 jag to 2D-IMAGE analysis (Figure 4). Although the 
diversity of molecules was shared among the different 
sources of DNA, some cells demonstrated a specific 
distribution of the major mtDNA topoisomers. For 
example, the slowest growing cell line MIN6, a pancreatic 
insulinoma line, harboured an additional band 
(Figure 4A, upper right panel, indicated by black dia- 
mond) and altered abundance of both ss and dsDNA 
structures, including low level of species 23 and 24. 

To determine how 2D-IMAGE profiles varied among 
tissues, which have a reduced mtDNA replication demand 
compared with actively growing cell lines, we isolated total 
DNA from four mouse tissues and two human tissues 
(Table 1). From mice, we analysed total DNA isolated 
from liver, heart, kidney and quadriceps (Figure 5A). 
As with cells, we observed the same major mtDNA topo- 
isomer species, but again there were clear differences 
in population distribution (Figure 5C). These differences 
are emphasized when total mtDNA copy number 
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Figure 4. 2D-IMAGE profiling reveals similar mtDNA structural forms with different distributions among many cell types. (A) Mouse cells: MEF, 
MIN6, C2C12 and MC3T3-E1. (B) Human cells: HMEC 151p (hTERT), MDA-MB-231, HeLa and SKBR3. MtDNA detected by Southern blot 
with total mtDNA probe. Black diamond in MIN6 cells indicates unique structure similar to those found in tissues (Figure 5). Abbreviations for 
structures are from Figure 3. 



(Figure 5D) is factored into their relative abundance 
(Figure 5E). In the two human tissues analysed 
(Figure 5B), pancreatic islets and left ventricle myo- 
cardium, we saw a clear laddering of catenated structures. 
These patterns are similar to results from others that 
focused on high molecular weight structures in human 
heart, such as junctional complexes and catenanes (27). 
2D-IMAGE analysis reiterated the high levels of caten- 
ated mtDNA molecules in not only human heart but 
also pancreatic islets. This heavily catenated mtDNA 
network has been suggested to be the replicating popula- 
tion in heart (28). As an extension of those findings, the 
simpler catenanes present in the majority of cell and tissue 
types shown here could be involved in replication or 
be replication products. Myocardium samples also 
contained the novel structure observed in MIN6 cells 
(black diamonds) (Figure 5B), which was also observed 
in approximately half of the mouse heart samples 
(Figure 5A). The relative distribution of the major 
mtDNA topoisomers differed between the two human 
tissues. Notably, human and mouse tissues showed a 
limited abundance of structures that contain ssDNA 
(structure #9, 13, 14, 18, 19, 23, 24 and 25). 2D-IMAGE 
profiling clearly revealed tissue and cell type-specific 
distributions of mtDNA topoisomers, suggesting that 
alterations in these profiles could be used to detect 
qualitative changes in mtDNA metabolism. 

Fully hybridized and partially hybridized RNA associated 
with discrete mtDNA structures in 2D-IMAGE 

RNA has long been known to be associated with the 
mtDNA nucleoid, playing roles in linking mtDNA to 
the membrane (41,8,7), DNA replication priming 



(42-44) and mtDNA conformation (37,45). To test 
whether RNA might be associated with specific DNA 
structures, we examined the effects of several 
well-characterized RNAses on the 2D-IMAGE profile. 
RNase Tl, which cleaves after unpaired G residues in 
RNA, and RNase VI, which cleaves base-paired RNA, 
had no overt effect on the 2D-IMAGE profiles relative 
to the untreated sample (Supplementary Figure S3). We 
were, however, able to detect changes in 2D-IMAGE 
profiles induced by treatment with RNase A, which 
digests single and double-stranded RNA after U and 
C bases in low salt (46), and with RNase H, which 
cleaves the RNA portion of a fully hybridized RNA- 
DNA duplex (Figure 6: dashed and solid arrows indicate 
decreased or increased signal after treatment, respect- 
ively). Relative to RNase A treatment (Figure 6C), un- 
treated samples showed significantly stronger 
hybridization, indicating the presence of more nucleic 
acid and with new hybridizing populations (Figure 6A: 
labelled a-e). Treatment with RNase H specifically 
removed one of these populations [Figure 6B (#), dashed 
arrow] to produce several additional populations (corres- 
ponding to Figure 3, #23-25). These RNase H-sensitive 
species may represent replication intermediates and were 
characterized further in subsequent experiments (see fol- 
lowing section). These effects on the distribution of signal 
were also shared by RNase A treatment (Figure 6C). 
Those structures generated by both RNase A and H 
were removed by SI nuclease treatment in sequential di- 
gestion experiments, indicating that they contain ssDNA 
(Figure 2F and Supplementary Figure S4). Our data show 
that RNase H-sensitive species are a subset of RNAse A 
(from several manufacturers) sensitive RNAs, with those 
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Figure 5. Mouse and human tissues share many common structures, but differ in relative abundance and additional structures. (A) 2D-IMAGE 
analysis of mouse tissues: liver, heart, kidney and quadriceps. (B) 2D-IMAGE analysis of human tissues: pancreatic islets and non-failed 
myocardium. Total DNA from indicated tissue was resolved by 2D-IMAGE and detected by Southern blot with total mouse mtDNA probe 
(panel A) or mixed human mtDNA probe (panel B). Black diamond indicates structure similar to that found in MIN6 cells in Figure 8. 
(C) Relative distribution of major topoisomers in four mouse tissues shown with standard deviation of biological replicates (± standard error of 
the mean) (n = 3). Asterisk indicates P<0.05; double asterisk indicates J P<0.01; triple asterisk indicates P< 0.001. (D) MtDNA copy number in 
mouse liver, heart, kidney and quadriceps as determined by Southern blot. Data from biological replicates (n = 4) are shown with standard deviation 
(± standard error of the mean). Liver mtDNA levels were used for normalization. Asterisk indicates same P-values as in panel (C). (E) MtDNA 
copy number-adjusted distribution of major topoisomers in tissues expressed in arbitrary units (A.U.). 



RNA species sensitive only to RNase A likely being par- 
tially hybridized to DNA. 

In untreated samples, signal surrounding the 2n caten- 
anes was much less distinct than in RNase A-treated 
samples [Figure 6 A (&)], with additional horizontal exten- 
sions observed to the right of In NC (c) and at the lower left 



of In NC (d). Hybridization within (e) and to the left of 
the CCC molecules (J) was significantly higher without 
RNase A treatment, suggesting that denaturation 
before transfer was not sufficient to remove all RNAse 
A-sensitive RNA. This qualitative pattern was seen in 
multiple independent preparations (data not shown). The 



e58 Nucleic Acids Research, 2013, Vol. 41, No. 4 



Page 8 of 14 




% y 



23* 4 24 25 ' 



RNase A 



/ 



D 7 

6 
5 
4 

® 3 



CO 



2 - 



-4 
-5 
-6 
-7 
-8 
-9 
-10 J 



o 
o 
o 



14 



U- -3 - 3-8 



J1 ' 



f 11 



15 



ET 



21 



rr 



, rL , n 



22 




■ RNase A 
□ RNase H 



Figure 6. Mitochondrial DNA topoisomers are associated with RNA in 2D-IMAGE profiles. (A) 2D-IMAGE profile of wild-type MEF DNA 
without RNAse treatment. New topoisomers relative to those in Figure 3 are assigned letters. (B) RNase H treatment, which digests RNA:DNA 
hybrids, reveals the vertical spikes of DNA that are sensitive to SI nuclease (23-25). Molecules that decrease on digestion are indicated with dashed 
arrows, whereas those that increase are indicated with solid arrows. Numbers correspond to topoisomers in Figure 3. (C) RNase A treatment, which 
digests heterogeneous RNA after U and C bases, reveals the typical 2D-IMAGE pattern shown in Figure 2. Dashed arrows indicate a reduction in 
signal, and solid arrows indicate increased signal. (D) Quantitation of change in abundance of topoisomers shown in panels B and C relative to 
untreated. 



signal within the diagonal region between In L and In NC 
species is strongly enhanced by the RNase A treatment 
(Figure 6C: #14). Assignments and signal intensity are 
plotted as fold change in Figure 6D. The majority of 
RNAs partially hybridized to DNA are associated with 
CCC mtDNA, which is the preferred substrate for tran- 
scription in bacteria (47). These data suggest that CCC 
molecules are a potential substrate for mitochondrial 
transcription. 

Discrete topoisomers contain heavy strand-specific ssDNA 

RNase H can remove RNA primers from replicating 
DNAs (48), implicating the RNase H-dependent 
structures in replication (49). Although these and other 
structures were sensitive to SI nuclease treatment, such 
sensitivity could be because of hairpin structures, nicks 
or abasic sites. To determine whether there were extensive 
single-stranded regions in the SI nuclease-sensitive 



structures, we examined the 2D-IMAGE hybridization 
pattern of non-denatured DNA (Figure 7). Alkaline de- 
naturation before transfer is required for hybridization of 
double-stranded nucleic acids with a single-stranded 
radiolabeled probe. Omission of this step has been used 
to characterize ssDNA-dependent replication of hyper- 
suppressive petite genomes in yeast (50). To this end, 
we resolved RNase A-treated DNA isolated from 
mitochondria by 2D-IMAGE in duplicate and transferred 
samples with or without alkaline denaturation (Figure 7A 
and B). RNAse A is required to eliminate RNA com- 
pletely from the profile to reveal the ssDNA species 
being tested (see Figure 6). In this system, the only detect- 
able mtDNA structures after non-denaturing transfer 
were the ln-2n molecules (Figure 3: #9), the lower 
termini of the vertical spikes of ssDNA emanating from 
In NC and In L molecules (Figure 3: #23, 24), a small 
portion of the lower part of the linear ssDNA vertical 
spike (Figure 3: #19) and less than unit-length 
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Figure 7. Strand-specific single-stranded mtDNA species resolved by 
2D-IMAGE. (A) Standard 2D-IMAGE profile of total wild-type 
(WT) MEF mtDNA transferred under denaturing (alkaline) conditions. 
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indicated molecules adhere to the positively charged membrane 
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mtDNA transferred under denaturing conditions and probed with 
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single-stranded linear DNA (<ln ssL; Figure 3: #25). 
Taken together with the SI nuclease sensitivity results, 
we conclude that these structures harbour significant 
ssDNA content, with an increase in dsDNA content as 
18 and 19 approach In NC and In L. Notably, these 
structures are most abundant in fast growing cell types 
in which mtDNA is rapidly duplicating. 

To determine whether these ssDNAs comprised a 
specific strand of mtDNA, RNase A-treated DNA 
isolated from mitochondria was resolved by 
2D-IMAGE, transferred with the denaturation step and 
probed with three sets of strand-specific oligonucleotides 
corresponding to different regions of the genome 
(Figure 7C-H and Supplementary Figure SI). The inclu- 
sion of RNase A ensured the removal of RNA and elec- 
trophoretically focused ssDNA into structures 18, 19, 23 
and 24. As expected, the catenanes, In NC, ln-2n, In L, 
and CCC molecules were detected by both heavy 
strand-specific and light strand-specific probe sets, as 
these structures are composed of dsDNA (Figure 7C- 
H). The heavy-strand probe sets (Figure 7C, E and G), 
but not the light strand sets (Figure 7D, F and H), 
detected In ssC (#23) and In ssL (#24), thus indicating 
that these structures comprise heavy-strand sequence. 
Light strand probes sets detected 18 and 19 near In NC 
and In L, confirming that these structures become pro- 
gressively more double stranded as they approach their 
fully double-stranded counterparts. 

Recovery from EtBr replication arrest identifies 
replication templates, intermediates and products 

Intact mtDNAs that contain extensive ss regions could 
be replication intermediates. To test this possibility, 
we perturbed mtDNA replication by EtBr arrest and 
release procedures. EtBr intercalates between base-paired 
nucleotides to disrupt the processivity of mitochondrial 
DNA polymerase, without significantly affecting nuclear 
DNA replication (51-53). We noted that 2D-IMAGE 
profiles of mtDNA isolated from wild-type MEF cells 
treated for 24 h with 100, 250 or 500ng/ml of EtBr 
showed a significant decrease in total hybridization 
signal in response to increasing EtBr concentration 
(Supplementary Figure S5A-D). Cells exposed to 500 ng/ 
ml EtBr contained <10% of the mtDNA signal present in 
untreated cells by qPCR or total counts in 2D-IMAGE 
(Supplementary Figure S5E), with all mtDNA structures 
affected (Supplementary Figure S5F and G). 

To restart mtDNA replication, we first exposed MEFs 
to 500ng/ml of EtBr for 24 h (no recovery; Figure 8A), 
and then released cells from arrest by washing out the 
EtBr with fresh medium. Cells recovered for 1, 3 or 6 h 
(Figure 8B-D) and mtDNA content were monitored by 
qPCR or total hybridization counts (Figure 8E). The 
2D-IMAGE profile of the 6-hour time point is similar to 
that of unexposed cells actively growing in culture. At this 
time point, elevation in mtDNA content was statistically 
significant (Figure 8E). Similar timing of the mtDNA 
increase was observed by qPCR in an independent 
experiment (Supplementary Figure S6). The majority of 
the radioactive signal arose from In NC molecules 
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Figure 8. Recovery from EtBr treatment stimulates replication of high-molecular weight mtDNA structures. (A) 2D-IMAGE analysis of wild-type 
MEFs cultured for 24 h in 500ng/ml of EtBr. Medium was changed to omit EtBr and cells harvested at 1 (B), 3 (C) or 6h (D) post-washout. 
(E) mtDNA content during washout as determined by qPCR (left) or total mtDNA hybridization signal for panels (A-D). (F) Distribution of 
signal across nine selected topoisomers as mtDNA levels recover. (G) Relative abundance of selected topoisomers from (F) after setting unrecovered 
sample equal to 1. (H) Model for template-product relationship for early and late recovery from EtBr-mediated mtDNA replication arrest. Solid 
arrows indicate template-product synthesis, and dashed arrows indicate conversion of one topoisomer to another. Location of nick in In NC is 
arbitrary. 
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(Figure 8F). We assessed the order of events during EtBr 
recovery by establishing changes in the relative abundance 
of various intact mtDNA structures (Figure 8G). At the 
early time point (1 h), increases in hybridization signal in- 
tensity were observed in high molecular weight mtDNA, 
such as catenanes (2n+ and 2n) and molecules that 'drag' 
out of the well (Figure 8B). These 'drag' molecules are not 
commonly detected in normal cells, and they only appear 
with 250 and 500 ng/ml of EtBr treatment (Supplementary 
Figure S5C and D); therefore, these molecules likely rep- 
resent highly interconnected tangles that cannot migrate 
into the gel efficiently because of their high molecular 
weight. In tightly correlated changes (Pearson's correl- 
ation: R = 0.99, ,P<0.01), 2n+ and drag molecules 
increased for the first 3h of recovery, but dropped to 
pre-exposure levels at 6h. These alterations are clearly 
visible in dark exposures (Supplementary Figure S7). 
Because replication templates are incorporated into 
product by semi-conservative replication, the depletion 
of these high molecular weight structures implicates 
them as the substrate for replication early in EtBr 
recovery. Also increasing after washout are 2n and 
scCCC, which increase 4-5-fold by 6 h washout, suggest- 
ing that these structures are replication products during 
the early stages of EtBr recovery. The increase in In NC 
was positively correlated with ln-2n abundance 
(Pearson's correlation: 7^=0.99, P<0.05), and at 6h, 
NCs were present at significant levels. Similar trends 
were observed in ssDNA molecule #18 at 6h. In 
contrast, In L and #19, both linear molecules, showed 
little if any relative increase at the 6-h recovery time point. 

Based on the timing and correlation of changes in 
the 2D-IMAGE profile, we propose a template-product 
relationship among several identified structures 
(Figure 8H). High molecular weight mtDNAs that accu- 
mulate during EtBr exposure are the primary template 
for replication early in the recovery and produce 2n cat- 
enanes and scCCC molecules, which could arise from 
dedicated replication pathways or simply as the result of 
differences in topology at completion of replication. 
Because ln-2n molecules are replication intermediates 
that originate from the In NC structure, these data 
suggest that late in recovery when these molecules 
increase-and likely during normal replication- In NC 
served as the template for replication. To reach this 
point, the early replication products (2n catenanes and 
scCCC) were likely converted to In NC species. 
Continuing replication from In NC, the complete replica- 
tion of both strands produced catenanes and scCCC, 
whereas a delay in initiation of second-strand synthesis 
would allow completion of first strand synthesis and dis- 
placement of ssC containing heavy-strand sequence 
(molecule #23). Initiation of second-strand synthesis on 
ssC formed a gapped circle (molecule #18) becoming pro- 
gressively more double-stranded during elongation to 
produce In NC. 

DISCUSSION 

Our 2D-IMAGE technique was built on the successes of 
ID separation of mtDNA structures to reveal multiple 



new forms of double-stranded circular mtDNA, including 
a range of supercoiled molecules, nicked circles and 
multiple catenated species, most of which are present in 
a variety of mammalian cells and tissues. This technique 
has also identified several forms of mtDNA that contain 
ssDNA structures, including ln-2n, In ssC and In ssL, 
which are likely replication intermediates. 2D-IMAGE 
differs from restriction digest-based 2D methods, such 
as 2D-AGE, which has been successfully used to 
examine replication and recombination intermediates 
(27,36,54,55). Although 2D-AGE uses digested DNA 
to detect low abundance replication or recombination 
structures, such as X and Y arcs, the topology of the 
intact genome containing these structures cannot be ascer- 
tained using that technique. 2D-IMAGE resolved ~25 
intact mtDNA species of diverse structural content, 
providing a survey of total mtDNA complexity. Even 
greater diversity in the mtDNA population can be 
expected by focusing on further resolution of specific 
structures [ex. (27)]. 

RNA in the nucleoid 

The assessment of RNA content in the nucleoid may elu- 
cidate mitochondrial transcription and genome stability. 
We found that most partially hybridized RNAs were 
associated with CCC mtDNA (structure e); given that 
highly supercoiled DNA plasmids are better transcrip- 
tional templates in bacteria (47) and adenovirus (56), 
our results are consistent with CCC mtDNA being a 
major template for transcription in vivo. This idea is sup- 
ported by observations of promoter-specific transcription 
on supercoiled templates using yeast mitochondrial RNA 
polymerase (POLRMT) in vitro (57), although mamma- 
lian POLRMT does not show promoter specificity (58). 
This contrast is not yet understood, and it suggests that 
additional factors may be contributing to transcription of 
mammalian CCC mtDNA in vivo. 

Alternatively, RNAs could play a functional role in the 
stability of mtDNA by linking the nucleoid to the mito- 
chondrial inner membrane (41). The mtDNA: membrane 
association (8,7) is likely important for coupled transcrip- 
tion-translation import (59), and by analogy to the 
bacterial nucleoid (60), has been suggested to be important 
in mtDNA replication (61), nucleoid division (62), segre- 
gation (63-65) and DNA repair (66-68). Furthermore, 
partially hybridized RNAs have been visualized by 
atomic force microscopy as bound to DNA, where they 
increase the apparent DNA writhe (45). In those DNA 
preparations, RNAs represent ~25% of the nucleotide 
content, have short 5'- and Spends and unmodified 
Spends distinct from total RNA transcripts that suggest 
a difference in processing and targeting relative to mature 
transcripts. Additional studies of the sequence and 
structure of these RNAs may yet reveal essential structural 
determinants of mtDNA nucleoid content and stability. 
Further supporting the idea that RNAs are likely part 
of the nucleoid is the presence of RNA-binding proteins 
in nucleoid preparations (69-72) and the requirement 
of some RNA-directed nucleoid proteins or their 
associated proteins, for mtDNA stability (73-75). 
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Replicating structures 

RNA priming is required for all current proposed models of 
mtDNA replication: coupled leading and lagging strand 
synthesis (54); strand-displacement or asynchronous repli- 
cation (76,77); and RNA incorporation throughout the 
lagging strand (44). The presence of discrete structures 
whose electrophoretic characteristics are altered by RNase 
H and/or SI nuclease treatment in 2D-IMAGE (structures 
18 and 19) suggests the presence of extensive regions of 
RNA:DNA hybrids. Similar stretches have been visualized 
on mtDNA by atomic force microscopy after treatment 
with RNase A and incubation with ssDNA-binding 
protein (55) and by analysis of replication intermediates in 
2D-NAGE (native agarose gel electrophoresis) gels (78). 

To determine how these structures respond to changes 
in mtDNA replication, we depleted mtDNA in MEFs by 
inducing replication arrest with EtBr, followed by release. 
Relative to human cell lines [e.g. 143B human sarcoma 
(54)], a much shorter time of exposure was required to 
deplete mtDNA levels (24 h), which coincided with a 
more rapid initiation of replicative recovery after EtBr 
washout. We observed high molecular weight mtDNA 
molecules accumulating under conditions of strong 
mtDNA depletion (500ng/ml of EtBr), which after 
release from arrest, began to increase in intensity first to 
produce 2n and scCCC molecules. At 6h post-washout, 
these high-molecular weight structures were diminished, 
coinciding with strong increases in In NC, ln-2n and 
structure 18. Based on these observations, we suggest 
that the early and late replicative mechanisms both 
produce 2n and scCCC molecules, but vary in template 
utilization. Specifically, early recovery may use the preva- 
lent >2n molecules as template, whose replication inter- 
mediates are not resolved in these gel conditions. Late 
recovery may use unit-length templates (In NC) to 
generate replication intermediates (ln-2n and structure 
18) and likely represents the mechanism occurring in un- 
challenged cells. 

New concepts generated by this study may prove 
important in refining our understanding of mtDNA repli- 
cation. For example, the significant differences in template 
utilization during EtBr-induced replication arrest and 
release have not been considered previously in the litera- 
ture. Moreover, the evident structural complexity of 
mtDNA is significant, and the separation of specific struc- 
tures will provide a clearer description of mtDNA features 
than will analysing the genome en masse. This latter bulk 
approach has been taken to describe strand-asynchronous 
and RNA incorporation throughout the lagging strands 
replication mechanisms (44), but additional refinement is 
required to investigate whether one or multiple mechan- 
isms exist specifically or simultaneously, and under what 
physiological conditions they are regulated. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Figures 1-7, Supplementary Methods, 
Supplementary Protocols and Supplementary References 
[79-81]. 
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